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Cooperative Spontaneous Emission as a Many Body Eigenvalue Problem 
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Institute for Quantum Studies and Department of Physics, Texas A&M Univ., Texas 77843 

(Dated: April 12, 2008) 

We study emission of a single photon from a spherically symmetric cloud of N atoms (one atom 
is excited, N-1 are in ground state) and present an exact analytical expression for eigenvalues and 
eigenstates of this many body problem. We found that some states decay much faster then the 
single-atom decay rate, while other states are trapped and undergo very slow decay. When size 
of the atomic cloud is small compared with the radiation wave length we found that the radiation 
frequency undergoes a large shift. 



Recent quantum optical experiments and calculations 
[U, Q focus on the problem in which a single photon is 
stored in a gas cloud and then retrieved at a later time. 
The directionality and spectral content of the coopera- 
tively reemitted photon is then of interest. 

Furthermore synchrotron radiation experiments in- 
volving N nuclei excited by weak 7 ray pulse have fea- 
tures in common with the present problem [S] . For exam- 
ple, in such experiments a thin disk of nuclei can easily 
be prepared in a superposition in which the atoms are 
all in the ground state together with a small probabil- 
ity of a uniform excitation of the state, similar to Eq. 
(|29)) . added in. The simplest example of two-atom coop- 
erative decay has been studied in many publications [J]. 
The A^-atom problem has been also investigated by sev- 
eral authors [5|] . Time evolution and directionality of the 
radiation emitted from a system of two- level atoms which 
are excited by a plane-wave pulse have been discussed in 

Having motivated our interest in the problem we now 
turn to the analysis of the correlated spontaneous emis- 
sion from N atoms in free-space. We consider a system of 
two level (a and b) atoms, initially one of them is in the 
excited state a and Ea — Eh = fko. Initially there are no 
photons. Atoms are located at positions Vj {j — 1, ..., A^). 
In the dipole approximation the interaction of atoms with 
photons is described by the Hamiltonian 
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_ (1) 
where dj is the lowering operator for atom j, a^ is the 

photon operator and gu is the atom-photon coupling con- 
stant for the k mode. We look for a solution of the 
Schrodinger equation for the atoms and the field as a 
superposition of Fock states 
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(3j{t)\bib2...aj...bN > |0> + 



number is equal to one and all atoms are in the ground 
state b. For simplicity we neglect the effects of photon 
polarization. Substitute of Eq. ^ into the Schrodinger 
equation yields the following equations for j3j (t) and 7^ (t) 
(we put h= 1) 

Pj{i) = -J^.9fe7fc(i)exp[-«(z/fe -uj)t + ik- r^], (3) 

k 



N 

7k(i) ^-i^gkPj{t)e-xp[i{vk-uj)t-ik-Vj]. (4) 
Integrating Eq. ([4]) over time gives [7/0 (0) — 0] 

r-t N 

7k(i) = -i / dt'^gkl3j{t')exp[i{vk-oj)t' -ik-Tj]. 
Jo .,=1 



Substituting this into ([3]) we obtain equation for /3j (t) 
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(6) 
We proceed by making the Markov approximation a-la 
Weisskopf and Wigner to obtain 



N 

where for i ^ j 

_ sin(fco|ri - Vj\) .cos(fco|ri - TjI) 



fco|r,-rj| 



fco|ri - r^l 



(7) 



(8) 



Tii = 1, fco = Lo/c and 7 is the single atom spontaneous 
decay rate 
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l\.{t)\bib2...bN > |lk > • 



(2) 



States in the first sum correspond to zero number of pho- 
tons, while in the second sum the photon occupation 



Vph is the photon volume. 

We point out that a rigorous treatment of the problem 
beyond the rotating wave approximation Hamiltonian P]) 
also yields Eqs. ^ and ([5]) [?, 's'] . Imaginary part of Fy 



in Eq. ([8]) appears due to a short range interaction be- 
tween atoms which is induced by electromagnetic field 
and causes a frequency shift Q. The frequency shift be- 
comes substantial when size of the atomic cloud is smaller 
then the wave length, this will be clear from Eq. ((26l) be- 
low. 

One can rewrite Eq. ([7]) in a matrix form 

B = -^TB, (9) 

where the vector B and the decay matrix T are given by 



B^ 



( Piit) \ 



r = r- = 



( 1 

r2i 



\PN{t) ) 



ri2 
1 



riAr\ 



The matrix T is symmetric Vi 



\ Tati Tn2 ■ ■ ■ 1 / 

(10) 
= Tji. Let |Ai > be 
eigenvectors of F and Xi {i = 1, ...,N) are the correspond- 
ing eigenvalues. If initially the system is prepared in an 
eigenstate -8(0) = jA^ >, then according to Eq. ^ the 
state evolution is given by 



B{t) = e-'-'^''B{0) 



(11) 



that is state evolves independently of other states and 
decays with the rate 7ReAi, where 7 is the spontaneous 
decay rate of a single atom. The general solution of the 
Schrodinger equation is 

* = Cie-^^i*|Ai > +C2e-^^^'|A2 >+...+ 
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\XN>+y27i.{t)\bib2...bN>\lk>, (12) 
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where Ci, C2, ..., Cn are constants determined by the 
initial conditions. Real part of Xi are positive numbers 
and, hence, the general solution corresponds to an expo- 
nential decay of the initial state. 

One should note that since tr(r) = N we obtain 
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^A.=tr(r) 
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i=l 



which is useful and insightful result. 

Next we solve the eigenvalue Eq. ([7]) analytically for 
a dense spherically symmetric cloud with atomic density 
p{r) {J p{r)dr = N). Tij changes in a scale of l//co- 
Assuming there are many atoms in the volume l/fep we 
can replace summation by integration, then Eq. ([7]) reads 
(/3 = e-T^*/3(r)) 



dr' p{r')[Ko{r,r') + iKi{r,r')] (3{r') ^ X(3{r), (14) 



where 
^o(r,r') 



sin(fcolr-r'l) 
fco|r-r'| 



Ki{r,r') 



cos(fco|r — r'l 
/colr-r'l 



It turns out that for koR ^ 1, where R is the char- 
acteristic size of atomic cloud, the eigenfunctions of the 
integral Eq. (fT4|l are determined by the real part of the 
kernel, Ko{r,r'), that is by equation 



dr'p(r')i^o(r,r')/3(r')=A/3(r). (15) 

Eq. ([15]) can be solved using the identity |10l | 
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and the orthogonality condition 



dnrY*„^{f)Yks{f) = SnkS„ 



where Ynm are spherical harmonics, f is a unit vector 
in the direction of r and j„(a;) are the spherical Bessel 
functions. Solutions of Eq. p5|l are 



Pnmir) = jn{kor)YnmiO, (fi) , 



Xr. 



drp{r)jlikor) 
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(17) 



where 9 and ip are angles describing direction of r in 
the spherical coordinate system. Eigenvalues (|17p are 
independent of to. That is the eigenvalues Xnm are {2n + 
l)-fold degenerate. If atoms are uniformly distributed 
inside a sphere of radius R (that is p{r) — N/V, for 
r < R) we obtain [l3| 

3iV 

Xnm = — [jnikoR] - Jn-l{koR)jn+l{koR)] • (18) 

In the limit koR '^ n we find 



/3„m(r) w - sin [kor - -nj Ynm{0, v), (19) 
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Fig. [T] shows Xnm/N as a function of n obtained from 
Eq. dUl) at different X/R < 1 (A = 27r/fco). The states 
with n < koR are degenerate and decay with the rate 
3jN/2{kQR)'^, while decay of the states with n > koR is 
suppressed. 

ii'i(r, r') term in Eq. (|14p yields an imaginary contri- 
bution in Xnm (frequency shift) and also modifies the real 
part. To solve Eq. ([T4|) we use the identity [iO| 

!5*^!^^4,if;i:n.(fw,(.'-') 

°' ' k=Os=-k 
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FIG. 1: Xnm/N as a function of n obtained from Eq. (|18|l at 
different X/R. 
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Jk{kor')h'^'{kor), r > r' 



jk{kor)h);>ikor'), r < r' 



(21) 
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where /ij. (z) are the spherical Bessel functions, h^, (z) = 

^7r/2ziJj,^^ ,2 (^)' ^k+i/2 (•^) ^^^ ^'^^ Hankel functions of 
the first kind. 

If atoms are uniformly distributed inside a sphere of 
radius R the answer for eigenfunctions is given by 



/3(i") = Jn {akor)Ynm{f), 



where 



a = Jl 
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0^ ^n 



(22) 



(23) 



and the eigenvalues A„ are determined from the following 
equation for a 



(1) 



jn{akoR) KUik^R) 



(24) 



jn-l{akoR) h'n\koR) ' 

In the Dicke limit fepi? < 1 Eqs. ^ and ^ yield 
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where A„; are nonnegative zeroes of the Bessel function 
jn-i {x). In particular, Aqi = {21 — l)7r/2 and An = ttI, 



/ = 1, 2, 3 . . .. One can see that Im(A„;) (frequency shift) 
becomes large for kgR -^ 0. In the Dicke limit only 
eigenvalues with n = have large real part and decay 
fast with the rate F; = 96N^/tt'^{21 — 1)^ (Dicke superra- 
diance [l2l | ) , while eigenvalues with n > are suppressed 
by a factor (fco-R)^". Those states are trapped. Please 
note that X^i^i ^i ~ -^7' ^^ expected from general ar- 
guments. To check validity of our analytical results we 
solved numerically the eigenvalue problem for the matrix 
(jlOp with atoms randomly distributed inside a sphere. 
In the Dicke limit our numerical simulations show excel- 
lent agreement with Eq. (pS)) already for a few hundred 
atoms. 

In the limit koR ^ n we find 
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(cos(2fcoi?) - 1) 



(27) 



In such a limit the contribution from the Ki term in the 
kernel is smaller by a factor 1/koR then those from the 
Kq piece. 

Figs. [2] and [3] show real and imaginary part of Aq/ as a 
function of koR obtained by solving Eq. (|24|) numerically. 
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FIG. 2: Real part of Aq; as a function of koR for I = 1, 2 and 
3. 

We next discuss the angular distribution of the emitted 
radiation produced by an eigenstate Pnim- Taking into 
account (3j{t) — e~^^"'-^(3nim{'>^j) we obtain from Eq. ([5]) 
for a uniform cloud of radius R 

7k(t) = 5k TT-^ — — : ^x 

I drr^ I dnrPnimir)exp[-ik-r]. (28) 




FIG. 3: Imaginary part of Xqi as a function of koR for I = 1, 
2 and 3. 



Substitute Eq. ^^ yields 7k (i) oc y„m(A:). That is the 
angular distribution of the emitted radiation is given by 
the same spherical function Ynm{0,ip) which describes 
anisotropy of the initial eigenstate. 

Let us next apply the proceeding to the problem of a 
uniformly excited sample in which one atom is excited 
and the state is given by 

1 ^ 
l*(0)) = -^E^'''"'^l^i'^2,...,a„...,M|0). (29) 

This is an interesting many body problem comprised of 
N atoms collectively emitting a single photon Jj. For 
example, this state can be prepared by the experiments 
of the type described in Ref. [4] . We are interested in the 
state of the system at time t in the limit koR ^ 1. We 
obtain the state evolution by expanding |^(0)) in terms 
of the eigenvectors | A„m > 

N 

|A„m >= ^jn{korj)Y,ijn{fj)\bi,b2, ■.■,aj, ...,6Ar)|0). 

3 = 1 

(30) 
We do this using the relation 



exp(iko-rj) = 47rE ^ i"j„(fcorj)y„*^(fco)>"nm(rj), 

n— m— — n 

(31) 
where fco and fj are unit vectors in the directions kp and 
Tj respectively. The state p9)) evolves as 
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where A„7 « 'i-iN/2{kQRf = Tjv, see Eq. ([2 
we have for the atomic state 



|*(i))atomic = e-^"*|^(0)) 



and therefore 
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g-TNtgiko-rj 
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Substituting this into ^ we obtain 



(34) 



-^ '^-^ - 3k [exp (-rjvt + i{yk - u;)t) - I] ^ „-»(ko-k).r, 
y/N to -Vk- »rjv ^ 

(35) 
which yields the following probability of a photon being 
emitted with the wave vector k 0, Q 



l7k(oo)| 
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The decay rate Fjv can be understood as the follow- 
ing. Dicke-like arguments [i2] for a coherent decay would 
yield the decay rate oiN^. However, the spontaneous de- 
cay rate is due to emission in all 47r sr directions, while 
the state ([29]) emits photon in a small diffraction angle 
\/R. As a result, the Dicke superradiance rate A^7 is 
reduced by the ratio of the solid diffraction angle }? jR? 
to 47r sr, this yields r^y. 

In summary, we studied correlated spontaneous emis- 
sion from N atoms in a spherically symmetric cloud 
which is a A^— body problem. For a dense atomic gas 
we obtained analytical expression for the eigenstates and 
eigenvalues of the system. We found that some states 
decay much faster then the single-atom decay rate, while 
other states are trapped and undergo very slow decay. 
In the limit X ^ R only states with n — decay fast 
(Fat ^ A7), all other states are trapped. In the op- 
posite limit A <C i? many states decay with the same 
rate Tn ~ N-fX^/R^. We also found that for A > i? 
the eigenvalues have a large imaginary part which corre- 
sponds to a frequency shift of emitted radiation. 
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